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1 INTRODUCTION 


ABSTRACT 


The level of star formation in elliptical galaxies is poorly constrained, due to difficul- 
ties in quantifying the contamination of flux-based estimates of star formation from unrelated 
phenomena, such as AGN and old stellar populations. We here utilise core-collapse super- 
novae (CCSNe) as unambiguous tracers of recent star formation in ellipticals within a cos- 
mic volume. We firstly isolate a sample of 421 z « 0.2, r « 21.8 mag CCSNe from the 
SDSS-II Supernova Survey. We then introduce a Bayesian method of identifying ellipticals 
via their colours and morphologies in a manner unbiased by redshift and yet consistent with 
manual classification from Galaxy Zoo 1. We find ~ 25% of z « 0.2 r < 20 mag galaxies 
in the Stripe 82 region are ellipticals (~ 28000 galaxies). In total, 36 CCSNe are found to 
reside in ellipticals. We demonstrate that such early-types contribute a non-negligible frac- 
tion of star formation to the present-day cosmic budget, at 11.2 + 3.1 (stat) 1 (sys) 96. 
Coupling this result with the galaxy stellar mass function of ellipticals, the mean specific 


star formation rate (SSFR; .S) of these systems is derived. The best-fit slope is given by 
log(S(M)/yr) = — (0.80 + 0.59) log(M/10!°5M.,) — 10.83 + 0.18. The mean SSFR 


for all log(M/Mo) > 10.0 ellipticals is found to be S=92+4+24 (stat) ee (sys) x107!? 


yr !, which is consistent with recent estimates via SED-fitting, and is 11.8 + 3.7 (stat) E 


(sys) 96 of the mean SSFR level on the main sequence as also derived from CCSNe. We find 
the median optical spectrum of elliptical CCSN hosts is statistically consistent with that of a 
control sample of ellipticals that do not host CCSNe, implying that these SN-derived results 
are well-representative of the total low-z elliptical population. 


Key words: galaxies: elliptical and lenticular, cD, transients: supernovae, galaxies: evolution, 
galaxies: luminosity function, mass function, galaxies: disc. 


ellipticals galaxies did indeed occur at high redshift. For instance, 
both rest-frame optical colours and the ‘Fundamental Plane’ of el- 


Elliptical galaxies dominate the high-mass end of the stellar mass 
function at late epochs (e.g. Kelvin et al. 2014), making them es- 
sential probes of our structure formation paradigm. The stellar mass 
growth of these systems is tightly linked to their star formation 
rates over cosmic time. The classical hypothesis of a ‘monolithic 
collapse’ (Eggen et al. 1962; Larson 1974; Chiosi & Carraro 2002) 
suggests that the stellar populations of early-type galaxies formed 
at high redshifts (z = 2), and that these galaxies have since evolved 
passively. 

Optical properties of ellipticals appear to obey simple scal- 
ing relations which indicate that the majority of star formation in 
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lipticals (Jorgensen et al. 1996; Saglia et al. 1997) exhibit only a 
small scatter, and do not appear to have evolved significantly for 
several Gyrs (e.g. Bower et al. 1992a; Franx 1993; Van Dokkum 
et al. 2000; Peebles 2002). 

Problematically, a monolithic star formation history cannot 
be aligned with the currently widely-accepted A-CDM framework, 
which in both semi-analytical models (e.g. Cole et al. 2000; Khoch- 
far & Burkert 2003) and in recent hydrodynamical simulations (e.g. 
Schaye et al. 2015; Pillepich et al. 2018) predicts ongoing contribu- 
tions to early-type galaxy stellar mass from merger and interaction 
driven star formation (Martin et al. 2017). Therefore, an accurate 
quantification of the star formation levels in ellipticals has ramifi- 
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cations for our knowledge of both early-type galaxy evolution and 
cosmology. 

The star-formation rates of early-type galaxies can be esti- 
mated observationally via several means. UV emission is a widely 
utilised indicator, as it exhibits the signal of photospheric light 
from the very youngest, most massive stars (e.g. Salim et al. 2007; 
Bouwens et al. 2009). A less direct indicator is mid and far IR 
emission, which captures the heating of dust due to young stellar 
populations (e.g. Calzetti et al. 2000; Inoue et al. 2000). Specific 
nebular recombination lines are also widely utilised, most notably, 
Ha, which has contributions to its flux from hot, luminous, young 
stars (e.g. Glazebrook et al. 1999; Erb et al. 2006). 

However, each of these diagnostics can have their signal con- 
taminated by other mechanisms. For example, post-AGB stars see 
their outer envelopes, in the process of being ejected, illuminated 
by core UV emission (Greggio & Renzini 1990). In Horizontal 
Branch stars, the intense temperatures of core-helium burning can 
also cause UV emission in these stars’ optically thin regions. In the 
case of both phenomena, strong emission lines such as Ha can also 
be generated; lines which are typically associated with the pres- 
ence of young stars (Preston et al. 2006). Other mechanisms which 
could be misinterpreted as star formation without careful investi- 
gation include the ionisation of gas due to active galactic nuclei 
(AGN; Groves et al. 2006), or due to shocks, with the collision of 
gas clouds (Oparin & Moiseev 2018). 

In the case of galaxies residing on the star forming main 
sequence (Noeske et al. 2007), such as most spiral galaxies, the 
contributions from these contaminants make up a small fraction of 
the star formation signal, and can be largely ignored. However, for 
galaxies below the main sequence, such as ellipticals, these con- 
taminants present a serious boundary to studies of their star forma- 
tion rates. 

A consequence of this issue is some degree of debate over the 
precise level of star formation typical in ellipticals, which in turn 
perpetuates the debate over our pictures of early-type galaxy mass 
assembly and cosmology. Whilst some studies conclude practically 
zero levels, (e.g. Jura 1977; Kennicutt 1998), 

more recent work suggests that low-level star formation per- 
sists in ellipticals over at least the latter half of cosmic time. For 
instance, Kaviraj et al. (2007) find that, while they are a homoge- 
neous population in optical colours (Bower et al. 1992b), ellipti- 
cal galaxies show more than 6 mags of colour spread in the NUV, 
which is a much wider range than can be produced from old stel- 
lar populations (e.g. Yi 2003). They conclude that at least 30% of 
local elliptical galaxies have NUV colours indicative of recent star 
formation within the last Gyr, and at z x 0.1, elliptical galaxies 
have 1% — 3% of their stellar mass in stars produced within such a 
time-frame. 

Interestingly, for ellipticals at intermediate redshift (0.5 < 
z < 1), ie. epochs at which the Universe is effectively too young for 
old stellar populations to exist, the characteristics of NUV colours 
remain the same (Kaviraj et al. 2008; Martin et al. 2018). This im- 
plies that this observed low-level star formation persists in early- 
type galaxies over at least the latter half of cosmic time. There also 
exists a strong correspondence between the presence of disturbed 
morphologies and blue UV colours (Kaviraj et al. 2011) indicat- 
ing that the star formation is merger driven (Schweizer et al. 1990; 
Schweizer & Seitzer 1992). Furthermore, the frequency of these 
morphological disturbances is larger than can be accounted for by 
the major-merger rate alone (Kaviraj et al. 2011), suggesting that 
minor mergers (i.e. mergers between ellipticals and gas-rich satel- 
lites) are likely to drive much of the elliptical star formation seen in 


the NUV, at least at late epochs (Kaviraj et al. 2009), with ~ 14% 
of the total cosmic star formation budget residing in the elliptical 
population in the low-redshift Universe (Kaviraj 2014a,b). 

Broad-wavelength SED-based derivations, such as those 
from Salim et al. (2016), offer arguably the most sensitive treatment 
to sub-main sequence star formation rates to date, due to the abil- 
ity to better identify sources of contamination to the star formation 
signal from the broad view of spectral information spanning from 
the UV to the IR. However, even these estimates are only reported 
as upper limits on the level of star formation in ellipticals. 

It is common to try to distinguish signatures of true star for- 
mation from those of AGN or old stellar populations using galaxy 
emission line ratios, such as Ho/[NII] and Hf/[OIII] (e.g. Kauff- 
mann et al. 2003). However, there is scatter in these relations, and if 
ellipticals do exhibit some level of star formation, they are unlikely 
to occupy the main locus as blue star-forming galaxies. As such, 
the issue is turned on its head, and low-level star formation could 
easily be misattributed to AGN activity and/or old stellar emission 
in a discrete treatment, such as that seen with the BPT diagram 
(e.g. Baldwin et al. 1981; Veilleux & Osterbrock 1987). Addition- 
ally, there are not an abundance of models which can predict these 
ratios from first principles. 

Recent sophisticated simulation suites such as EAGLE 
(Schaye et al. 2015), Horizon-AGN (Dubois et al. 2016; Kaviraj 
et al. 2017) and IllustrisTNG (Pillepich et al. 2018) support the 
suggestion that minor mergers could generate a persistent low-level 
of star formation. This would not only make them a phenomenon 
crucial to the evolution of ellipticals, but would also point further 
towards the hierarchical evolution now widely accepted under a A- 
CDM paradigm. 

Such simulations also implement sub-grid star formation 
recipes as well as feedback mechanisms which can remove cold 
gas and quench star formation. Since these sub-grid recipes are 
often calibrated and benchmarked against a number of observed 
properties and relations (including star formation rate densities; e.g. 
Pillepich et al. 2018), any biases or systematic uncertainties in these 
quantities can produce unrealistic sub-grid physics and degrade the 
predictive power of the simulations. 

If the observed elliptical morphology and colour are repro- 
duced organically via ‘dry’ mergers as a consequence of quenching 
through feedback, both merger rates (minor and major) and feed- 
back mechanisms could be better constrained in light of accurate 
and unbiased tracers of elliptical star formation rates. Furthermore, 
since many of these simulations now aim to replicate the star for- 
mation rate density (SFRD) of main-sequence galaxies as standard, 
a natural progression would be to test model physics using the ob- 
served star formation rates of the early-type population. 

There are other mechanisms which may also drive star forma- 
tion in ellipticals. For example, galaxies of all morphological types 
are surrounded by significant reservoirs of cold gas (Chen et al. 
2010; Thom et al. 2012), which may fuel low levels of star forma- 
tion even in apparently dead ellipticals (Tumlinson et al. 2017). In 
hydrodynamical simulations, the structure and physical properties 
of the CGM are found to depend strongly on factors such as res- 
olution, feedback and the presence of magnetic fields and cosmic 
rays (Hummels et al. 2019; Butsky et al. 2020). Therefore, accurate 
measurements of the level of star formation in elliptical galaxies are 
an important clue for understanding the connection between galaxy 
evolution and the CGM. In more sophisticated simulations, which 
incorporate magnetic fields and cosmic ray heating in addition to 
stellar and AGN feedback processes, observations like these will 
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be important for placing constraints on the recipes used to model 
these processes. 

For each of the described reasons, we are encouraged to 
search for an independent diagnostic of star formation which is not 
subject to the aforementioned uncertainties. Core-collapse super- 
novae (CCSNe) are unambiguous indicators of recent star forma- 
tion. The short (~ 6-40 Myr; Smartt 2009; Botticella et al. 2017) 
lifetimes of their progenitor stars, relative to the timescales of 
galaxy evolution, mean that the presence of CCSNe in the low-z 
elliptical population would offer explicit evidence of recent low- 
level star formation within them. 

These transients also provide a more thorough analysis of star 
formation in the shape of the volumetric SFRD (psrg) as a function 
of galaxy stellar mass, M. It was shown in Sedgwick et al. (20192) 
and Sedgwick et al. (2019b) that the SFRD of low-surface bright- 
ness galaxies can be accurately traced from the volume's CCSN 
rates, Occsn, as measured by the untargeted, high-cadence, SDSS- 
II Supernova Survey (Sako et al. 2018), by matching CCSNe to 
their host galaxies. This connection is shown via Equation 1, where 
f is the mean ratio of CCSN events to mass of stars formed. 


PsrROM) = peess (V) 


(1) 

In the case of low surface brightness galaxies (LSBGs), CC- 
SNe not only helped locate galaxies which would have been missed 
using traditional galaxy extraction pipelines, by acting as signposts 
towards them, but they also allowed the population's star forma- 
tion to be quantified in situations where spectral and broad-band 
photometry-based derivations would be subject to large uncertain- 
ties related to low signal-to-noise. It is this latter point which partic- 
ularly encourages us to extend a similar method to elliptical galax- 
les. 

The SFRD would allow us to place constraints on the ellipti- 
cal contribution to the cosmic star formation budget. As well as the 
SFRD, we may also be able to assess the mean specific star forma- 
tion rate (SSFR; S) of ellipticals, both as a function of mass, and 
for the total population. This is due to the connection between the 
SFRD, the galaxy stellar mass function (GSMF) and (S), demon- 
strated via Equation 2, where ¢ is the number density of galaxies in 
a cosmological volume. 


PsFr(M) 
QM) M 


In Sedgwick et al. (2019b), the SFRD was first found from 
CCSNe, before the GSMF was determined using a re-arrangement 
of Equation 2 from the assumption that S is constant with mass 
for low-surface brightness galaxies, due to their residence on the 
star forming main sequence (McGaugh et al. 2017). In the case of 
massive ellipticals, however, the GSMF is already well-constrained 
via various experiments (Franceschini et al. 2006; Vulcani et al. 
2011; Kelvin et al. 2014), and it is instead S which is our unknown. 
Whilst we will here utilise Equation 2 to measure the mean SSFR 
of ellipticals, the relation can be applied to any well-defined galaxy 
population for which any two of pspr, ¢ and SS are known. 

The structure of the present work is summarised as follows: 
Sections 2 and 3 summarise the definition of the CCSN and ellip- 
tical samples, respectively along with summaries of the data sets 
used. Section 4 gives measurements of the volumetric star forma- 
tion rate density of ellipticals. Section 5 shows results for the galaxy 
stellar mass function and mean specific star formation rate of ellip- 
ticals. Section 6 compares the spectra of CCSN-hosting ellipticals 


S(M) = (2) 
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to that of a control sample of standard ellipticals. Finally, Section 7 
summarises each of our findings in the context of early-type galaxy 
evolution. 


2 THE CCSN SAMPLE 


We begin with the SN sample utilised in Sedgwick et al. (2019b). 
This was formed from the SDSS-II Supernovae Sample (Sako et al. 
2018), cross matched with the IAC Stripe 82 legacy sample (Fliri & 
Trujillo 2016). A more detailed discussion of this cross-matching 
is given in Sedgwick et al. (2019b), which describes the careful 
assignment of host galaxies, and the rejection of variable stars and 
AGN. This sample was cut to exclude SN of r-band peak brightness 
» 21.8 magnitudes, fainter than which the sample is estimated to be 
incomplete (as also described in the aforementioned study). This 
initial sample consists of 2528 rj4;; < 21.8 mag SNe and their 
hosts, following variable star/AGN removal. 

We utilise the best redshift estimate available for each 
SN/host pair: this is spectroscopic, from either SN (taking pref- 
erence) or host galaxy, the latter from either SDSS-II Legacy (York 
et al. 2000), SDSS-II Southern (Baldry et al. 2005), SDSS-III 
BOSS (Dawson et al. 2013) or SDSS-IV eBOSS (Dawson et al. 
2016), or in the absence of spectra, photometric, from ugriz Stripe 
82 legacy survey galaxy fluxes, using the ‘scaled flux matching’ 
technique (Sedgwick et al. 2019a; Baldry et al. 2021). The SN sam- 
ple is then cut to z « 0.2, leaving 1070 SNe, 845 of which have a 
spectroscopic redshift. 

In the present work, we treat SN classification differently 
than in Sedgwick et al. (20192) and Sedgwick et al. (2019b). Sako 
et al. (2018) attempted to classify each of these SNe as either a 
Type Ia, Type Ib/c, or Type II. As a brief summary, 2 main sets of 
probabilities were calculated for each SN: 

The first was a set of Bayesian probabilities of belonging 
to each of the 3 SN classes, obtained by analysing observed light 
curves with the Photometric SN IDentification (PSNID) software 
(Sako et al. 2011), using a grid of Type Ia, Type Ib/c and Type II 
SN light curve templates. 

The second was a set of Bayesian + nearest-neighbour (NN) 
probabilities, the improvement being that with the NN extension, 
probabilities account for differences in the distributions of extinc- 
tion, light curve shape and redshift seen for each of the 3 SN types. 
Synthetic SNe of each type were considered a ‘near-neighbour’ 
to an observed SN if within a threshold Cartesian distance in a 
3-dimensional parameter space consisting of the above 3 param- 
eters. The NN probability of the candidate being a Type II SN, for 
instance, was then the fraction of NN simulated SNe which were 
Type II’s. In the present work, we use the NN probability, Pyyn if 
there are at least 20 near-neighbours for the transient, which is the 
case for 738 z « 0.2 SNe, and in the 332 remaining cases we use 
the Bayesian-only probability, Pgayes. We denote the best probabil- 
ities available as Pia, Proc and Py. Though there is some spread in 
the Pyyn — Ppayes distributions, with a standard deviation of 0.22 for 
the z « 0.2 Type Ia probabilities, for example, the mean difference 
of Pyn — Peayes for Type Ia SNe is very close to zero, implying a 
consistency between the 2 methods on a statistical level. 

Elliptical galaxies are more traditionally associated with 
Type Ia SNe than with CCSNe (e.g. Gallagher et al. 2008). To 
ensure Type Ia interlopers do not cause an overestimate of star- 
formation properties, we therefore require Pj, < 0.05 for the SN to 
be included in our CCSN sample. Removing those with Pia > 0.05, 
we are left with 360 confidently classified z < 0.2 CCSNe. This 
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sub-sample is used for all main CCSN-based results in the present 
work. However, we will repeat our analysis including an additional 
61 0.05 < Pia < 0.5 objects, to test for the sensitivity of results to 
Type Ia contamination. Note that these numbers suggest the CCSN 
probabilities are largely bi-modal. Put differently, ~ 93% of all 
z < 0.2 SNe have either Pj, < 0.05 or Pj, > 0.95. 

The ratio of likely Type Ib/c SNe to likely Type II SNe is also 
approximately in line with expectations from the literature (e.g. 
Hakobyan 2008), with 24 + 1 (Poisson error) % having Pr. > Pr. 
Furthermore, it is far more common to misclassify a Type Ia as a 
Type Ic rather than as a Type II SN, due to the similar light curve 
shapes of Type Ia and Type Ic supernovae (lack of plateau), and 
due to their similar spectra (Clocchiatti & Wheeler 1997). As a re- 
sult, if we had significant contamination from Type Ia supernovae 
in our core-collapse sample, we would expect a higher ratio of Type 
Ib/c to Type II supernovae than is observed. These points increase 
our confidence that the Bayesian SN-type classifications are fully 
trustworthy. 

Other tests for the sensitivity of forthcoming results to con- 
tamination include a repeating of results with SN type probabili- 
ties re-scaled to account for the efficiency and purity estimates in 
each SN class. These values were estimated from a set of simulated 
SNe by Sako et al. (2018). We also test the effect of simply us- 
ing Bayesian-only probabilities for all SNe. As might be expected 
given the aforementioned mean of the Pyy — Ppayes distribution, and 
the bi-modality of the CCSN probabilities, respectively, we find for 
both tests that our forthcoming results do not change notably. In 
conclusion, the CCSNe are classified from their light curves with 
sufficient confidence that our sample selection is stable to the fine 
tuning of the SN classification procedure. 


3 DEFINING ELLIPTICAL GALAXIES 


We must next determine which hosts galaxies are ellipticals, and as 
quantitatively as possible. For this we turn to the wider sample of 
galaxies in Stripe 82 legacy coadded imaging. 

Using estimates by Fliri & Trujillo (2016) from simulated 
de Vaucouleurs profiles, we can expect 95% of all z « 0.2 bulge- 
dominated galaxies in the legacy coadded imaging to have at least 
95% of their light recovered at r,,4; = 20 mag, where auto denotes 
the near-total light ‘Kron’ elliptical aperture (Kron 1980; Bertin & 
Arnouts 1996). We therefore impose this magnitude cut, leaving a 
sample of 113239 galaxies. Note that we do not apply this cut to 
the SN hosts, as it is the SN light which determines their detection. 
Nonetheless, 70% of all z « 0.2 SN hosts in our sample have r < 20 
mag, and this rises to 98% for M > 10'°°M., masses below which 
we expect few ellipticals. 

The most common and straightforward way to attempt to iso- 
late ellipticals is to use a cut on galaxy morphology. For exam- 
ple, Kelvin et al. (2014) define ellipticals as all single-component, 
bulge-dominated galaxies. A useful parameter to quantitatively rep- 
resent ‘bulge dominance’, available for all SDSS galaxies, is the 
‘de Vaucouleurs fraction’ (fjey). In SDSS, the composite model r- 
band flux for a galaxy, Fomp, is calculated by taking a best-fit linear 
combination of flux in a de Vaucouleurs profile with that in an ex- 
ponential profile, such that Feomp = faevFaev + (1 — faev)Fexp, Le. 
Feomp = Faev and fay = 1 for the most bulge dominated cases. 

Whilst defining ellipticals based on this parameter alone is 
sufficient for most applications, if we are to estimate the typical star 
formation of ellipticals we must take particular care to avoid any 


contamination from non-ellipticals, which could have a significant 
weighting on such calculations. 

To estimate the contamination when classifying on fav 
alone, we cross-match, within 2.5", the galaxy sample with Galaxy 
Zoo 1 (combined main and bias studies; Lintott et al. 2011), be- 
fore isolating those z « 0.2, r « 20 mag galaxies (although most 
have r < 18 mag) with Pgaz > 0.8 (4692 galaxies) and those 
z < 02, r < 20 mag galaxies with Pgz < 0.2 (3874 galaxies), 
where Pegaz is the fraction of votes for an elliptical classification 
within the Galaxy Zoo Project. 

Of those confidently classified galaxies with either Pgaz > 
0.8 or Pgaz < 0.2, we find that 20% of galaxies with faey > 0.5 
have Pggz « 0.2. Imagine this interloping 2046 possess (a maybe 
even conservative) 10 times the star formation rate of the genuine 
ellipticals in the sample: Suddenly, ~ 70 of all CCSNe (and star 
formation) within this sample would stem from interlopers. 

A naive solution would be to inspect the ~ 45000 z « 0.2, 
r < 20 mag, fav > 0.5 galaxies manually to remove non-ellipticals. 
A consistent inspection of this number of galaxies is not only im- 
practical but we also find that an over-reliance on manual inspec- 
tion results in a significant redshift bias towards the selection of 
nearby objects, due to the better resolution of galaxy structure at 
lower redshifts. 

We instead find that a Bayesian classification method, effec- 
tively based on distributions of colour and morphology, yields a 
near-complete sample of ellipticals, whilst removing contaminants 
effectively, and minimising redshift selection biases. 

We treat our previously obtained samples of Pg oz > 0.8 and 
Pggz < 0.2 galaxies as respective reference samples of ellipticals 
and non-ellipticals. We use redshift-debiased values of Pg cz from 
the Galaxy Zoo project's spectroscopic sample, and the raw clas- 
sification fractions from their photometric sample. Note, however, 
that for the low redshift range we are concerned with in the present 
work, the median change to classification fraction is small, at ~ 5% 
for 0.15 < z < 0.2, r < 20 mag galaxies. 

We next arrive at an elliptical classification confidence level, 
Pg Bayes; using a range of variables i, i + 1,i  2,..., n. Let P(E | i) 
correspond to the probability that a galaxy would belong to the 
Pggz > 0.8 sub-sample, given the galaxy's value for a variable i. 
Likewise, let P(N | i) denote the probability that the galaxy would 
belong to the Pgaz < 0.2 sub-sample given this input. According 
to Bayesian statistics, the odds ratio P(E | i)/ P(N | i) is then given 
by Equation 3, where P(E)/P(N) is the odds ratio prior to the con- 
sideration of our variable i. Crucially, P(i | E)/P(i | N) equates to 
the ratio of probability density function heights at the value for i. 


P(E|i  PG|E) P(E) 
P(N |i) PGI N) P(N) 

We assume as an initial prior that P(E)/P(N) = 1. As the 
odds ratio becomes the new prior when we combine information 
from successive variables, it follows that a galaxy’s elliptical clas- 
sification confidence level is given directly from the product of sev- 
eral PDF height ratios, as shown by Equation 4. 


(3) 


E -1 
pee or (4) 
P |N) 


i 


P E,Bayes = 


We use 10 variables in total to calculate Pg Bayes, as shown in 
Figure 1. The PDFs used as inputs for Equation 4 equate to Gaus- 
sian kernels fitted to each variable’s distribution, each smoothed 
empirically to avoid non-physical discontinuities. 
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Figure 1. The 10 sets of probability density functions used to classify elliptical galaxies (see text for details). Dark-red series denote the training sample of 
Pggz > 0.8 z < 0.2 r < 20 mag galaxies in the Stripe 82 region. Blue series show the same but with instead Pg az < 0.2. Gaussian KDEs are fitted to the 
distributions. For a given galaxy, the ratios of the KDE heights at a given value for each variable are used to calculate Pg Bayes. 


We firstly use rest-frame ugriz colours within elliptical Auro 
apertures, corrected for Galactic extinction using the maps of 
Schlegel et al. (1998) and k-corrected with the prescription of 
Chilingarian et al. (2010). These primarily help isolate 'red se- 
quence’ galaxies from the star-forming ‘blue cloud’ (see, e.g. 
Baldry et al. 2004). 

We next use the difference between auto aperture colours and 
those within 2.5" radius circular apertures (denoted APER, Bertin & 
Arnouts 1996). Late-type galaxies typically feature disks which are 
bluer than their bulges, due to the higher star formation levels in 
the former regions. Conversely, ellipticals display relatively radi- 
ally consistent colours. As a result, the wider spread about a null 
colour difference seen for the non-elliptical reference sample helps 
us to exclude disk galaxies. 

We finally use 2 measures of morphology. The first is the 
apparent ellipticity, measured in the r-band. The second is fiy. 
These 2 morphological parameters primarily aid the exclusion of 
dusty, red, star-formers. 

The distribution of final PẸ Bayes values found from the 10 
sets of input PDFs in unison is shown as the black series in Fig- 
ure 2. Approximately 95% of all z « 0.2, r « 20 mag galaxies have 
Pggayes > 0.95 or PE Bayes < 0.05. This includes all galaxies, not just 
the training sample. 

Also shown are the distributions of PẸ Bayes values resulting 
from the parameters used in isolation. It can be seen that galaxy 
colours are the most crucial parameters for a confident elliptical 
classification, with ~ 90% of z « 0.2, r < 20 mag galaxies having 
PẸ Bayes > 0.95 or PẸ Bayes < 0.05 from colour alone. Morphology 
(fav and ellipticity) and AuTO-APER colour differences are of sec- 
ondary and comparable importance for classification confidence. 

Of course, confidence is not always a reflection of the ac- 
curacy of classifications: Each of the parameters helps remove a 
different sort of contaminating object, and must be used in unison 
for an effective isolation of ellipticals. This is emphasised with the 
Stripe 82 coadded images included in Figure 2, which show galax- 
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ies that would have been assigned PẸ Bayes > 0.5 were it not for a 
given criterion. 

The top image shows a nearby irregular galaxy which ex- 
hibits similar morphological properties to a massive elliptical, but 
was rejected after a consideration of its optical colours. The central 
image shows an edge-on, dusty, star-forming galaxy which is con- 
sistent in colour with an elliptical, but was rejected due to faey and 
ellipticity. The bottom image shows either a poorly resolved spiral 
or a lenticular disk, consistent in both colour and morphology with 
an elliptical, yet rejected due to its AuTO-APER colour differences. 

Of a total of 113239 z < 0.2, r < 20 mag galaxies, 27940 
(^ 25%) have PE Bayes > 0.5, and it is these galaxies which define 
the elliptical sample, to be used in successive sections of the present 
work. 

Manual inspection of these objects reveals little to no con- 
tamination from any obvious non-ellipticals, nor any obvious evi- 
dence of LSBGs along the line-of-sight to these ellipticals which 
could have instead housed the observed CCSNe. Figure 3 shows 2 
more quantitative tests for the validity of our classifications. The 
left-hand panel shows k-corrected r-band absolute magnitude ver- 
sus rest-frame (g — r) colour, where all magnitudes are in auto 
apertures, with the aperture defined in the r-band. Contour levels 
are 10% to 90% of the peak number density in steps of 10%. The 
fact that the dashed-red and dashed-blue contours overlap shows 
that, much like a rudimentary cut on fgey, a ‘hard’ cut on colour 
would not be able to remove non-ellipticals without also removing 
a significant number of ellipticals. 

We see that our sample of ellipticals resulting from the 
Bayesian classification method (solid-red) follows a similar distri- 
bution to the training set of Pg.gz > 0.8 galaxies (dashed-dark-red), 
with the only notable difference being an offset in absolute magni- 
tude of ~ 1 mag between the peak density of the 2 distributions. 
This is due to the bias of having more confident manual classifi- 
cations in Galaxy Zoo for brighter galaxies. This is also clear from 
the fact that the magnitude at which we find peak number density in 
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Figure 2. Distributions of Pg Bayes values for the sample of ~ 10° z<02r< 20 mag Stripe 82 galaxies. In black: the final Pg paye, values using all 10 
variables in unison. In red, grey, and magenta: PẸ Bayes values from galaxy colours, morphology (ellipticity and fgey), and AuTO-APER colour differences alone, 
respectively. Adjacent images show examples of galaxies with final values of PE Bayes < 0.5, but which would have had PE Bayes > 0.5 were it not for a single 
set of variables, where image border colour denotes the set of variables which caused the crucial reduction to Pg Bayes. 
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Figure 3. Left: Rest-frame r-band absolute magnitude versus rest-frame (g — r) colour, where all magnitudes are those within elliptical auto (‘Kron’) apertures 
and are corrected for Galactic extinction. Contour levels are 10% to 90% of the peak number density in steps of 10%. Dashed-dark-red and dashed-blue 
contours depict the training samples of Pg az > 0.8 and Pg cz < 0.2 galaxies, respectively. The solid-red contours represent the sample of ellipticals resulting 
from our Bayesian classification method. Black points show the distribution of all z « 0.2, r « 20 mag Stripe 82 galaxies. Right: Galaxy counts as a function 
of specific-star formation rate as measured by GSWLC-2. The dashed-dark-red, dashed-blue, red-filled and black-dotted series represent the same respective 
sub-samples, described for the left-hand panel. The red-dotted series depicts the same sample as shown by the red-filled series, but excluding the training set 
of Ppz > 0.8 galaxies. 
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the Pg oz < 0.2 sample (dashed-blue) does not correspond the mag- 
nitude at which peak density is seen for all galaxies (black points). 

As a comparison with the CCSN-derived star formation re- 
sults of the present work we will repeatedly compare with star 
formation rates derived from UV/Optical SED fitting applied to 
the second edition of the GALEX-SDSS-WISE Legacy Catalogue 
(GSWLC-2; Salim et al. 2016, 2018). 

In the present work, 23009 out of 113239 z « 02 r « 20 
mag galaxies are matched with the GSWLC-2 sample within 2.5". 
These matches are ry; < 18 mag galaxies with a spectroscopic 
redshift in the range 0.03 « z « 0.2, which lie within the GALEX 
footprint (Martin et al. 2005; Morrissey et al. 2007). 

Star formation properties of these galaxies were estimated 
using SED fitting with the Code Investigating GALaxy Emission 
(cIGALE; Noll et al. 2009; Boquien et al. 2019). Excellent summaries 
of the SED fitting technique are given in Salim et al. (2018) and 
Turner et al. (2021), but in short: Synthetic spectra, generated us- 
ing the simple stellar population templates of Bruzual & Charlot 
(2003), based on a Chabrier (2003) initial mass function and cover- 
ing a wide range of metallicities (log(Z) = —2.4 to —1.3), are fitted 
to the observed UV-to-Optical photometry. 

Templates were then combined with Myr-resolution star for- 
mation histories. The library of these histories were built using ex- 
ponentials with various decay times for an old stellar population, 
with a relatively flat burst superimposed (at least 100 Myr ago) to 
represent a younger population. The SED estimates of star forma- 
tion rate were additionally constrained by the galaxy IR luminosity, 
itself calculated by matching the energy absorbed by a galaxy's dust 
with the energy it re-emits. Other properties derived in this fitting 
procedure include (but are not limited to) galaxy stellar masses, 
stellar ages, stellar metallicities, absolute magnitudes and colour 
excess. In the present work, we use the SFRs derived using the 
deepest photometry available for each galaxy (catalogue GSWLC- 
X2). 

The right-hand panel of Figure 3 shows the matched distri- 
bution of SED-derived SSFRs (.Ssgp). It is seen that the peak value 
of Ssgp comes at ~ 107129 yr! for both the training sample and 
Bayesian-classified sample of ellipticals. A notable difference be- 
tween our final elliptical sample and the Galaxy Zoo training sam- 
ple is that the former exhibits a secondary peak at Ssgp ~ 10-!!9 
yr !. However, we find no significant difference in the appearance 
of the galaxies classified as ellipticals at these 2 different levels of 
star formation, and so the origin of this secondary peak remains 
unexplained. 

According to our classifications, 74% of Ssep < 10-!!9 
yr are classified as ellipticals. This number rises to 80% for 
Ssep < 10-29 yr. Ellipticals become the dominant class for 
Sskp < 10-108 yr. 

In summary, this Bayesian method is able to isolate a near- 
complete sample of z « 0.2, r « 20 mag elliptical galaxies in the 
Stripe 82 region. In theory, this method could be extended to effi- 
ciently classify several species of galaxies over wide-field survey 
footprints, in a manner consistent over redshift. 


1 


4 THESTAR FORMATION RATE DENSITY IN 
ELLIPTICALS 


Our CCSN hosts were classified simultaneously with the larger 
galaxy sample, and as such are subject to identical classification cri- 
teria as in the previous section. Of our sample of 421 likely CCSN 
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hosts (Pia < 0.5), 36 are classified as ellipticals using our Bayesian 
procedure. 27 of these have Pra < 0.05. 

Co-added Stripe 82 images of these 36 galaxies are shown in 
Figure 4. These host galaxies can now be used to estimate the con- 
tribution to the cosmic star-formation density from ellipticals. We 
can utilise Equation 7 to estimate the CCSN-rate density per unit 
volume V in ellipticals using the sample of elliptical CCSN host 
galaxies (see also Sedgwick et al. 2019b). Galaxy stellar masses 
are first calculated using the same prescription as in Sedgwick et al. 
(2019b); a k-correction inclusive prescription based on i-band auto 
mag, (g—i) observed colour (in 2.5" circular apertures) and redshift. 

To correct for the fact that our SN sample is flux-limited 
but we instead want volume-limited CCSN statistics, we imple- 
ment a statistical correction identical to that used in Sedgwick 
et al. (2019b). This is similar to a 1/Vinax method, but where the 
SN light is the determining factor in detection, not the host galaxy 
light. We employ the volume-limited absolute magnitude distribu- 
tions of Richardson et al. (2014) for Type Ib/c and Type II SNe. 
For each SN, i, in our own sample, we use the mean of the abso- 
lute magnitude distribution for its SN type j, denoted M j. We can 
then compute the mean expected apparent magnitude for our SN, 
mi, as a function of its redshift (z;), Galactic extinction (A,w.;), 
host-galaxy extinction (A,4;) and k-correction (k,;), as shown by 
Equation 5. 


m; = Mj-*5logdi(zi) + ki + Arni + AcMwi (5) 


Using the standard deviation in the relevant absolute magni- 
tude distribution, c; we can estimate the detectable fraction, ej, of 
SNe that would have r « 21.8 mag, for the redshift, extinction and 
k-correction of our SN, using Equation 6. The inverse of the frac- 
tion ej; is the SN's (and hence its host's) weight of contribution to 
the volumetric number density. 


] 1 m; — 21.8 
ini ca Vio, | (6) 
B 1 >»; 1/&(M) 
PccsnM) AlogM 7V (7) 


Equation 7 shows that summing over a mass bin of width 
Alog M, the quantity >); 1/&(M) leads to the volume-corrected 
number of CCSNe associated with galaxies for each bin. t is the 
effective rest-frame time over which CCSNe could be identified by 
the survey. This is shorter than the observed time-frame of the un- 
targeted supernova survey, t, such that t = t / (1+z). The observable 
volume, V, is derived from the sky coverage of the Stripe 82 region 
and the imposed z < 0.2 limit. The relationship between the CCSN 
rate and SFRD is then given by Equation 1. As derived in Sedgwick 
et al. (2019b), we adopt logR = —1.9. 

We utilise a 10° iteration Monte Carlo (MC) technique to 
account for uncertainties in redshift, galaxy magnitudes and the 
resultant stellar masses which all feed into our calculation of the 
CCSN-rate density and SFRD as a function of galaxy stellar mass, 
the results for which are shown in Figure 5. The statistical errors 
(labelled ‘stat’) on these densities are represented by the filled re- 
gions, and equate to the quadrature sum of MC errors (standard 
deviation of the densities over MC iterations), Poisson errors and 
the cosmic variance, the latter estimated given the volume of the 
Stripe 82 region out to z = 0.2 (Driver & Robotham 2010). The red 
series denotes the estimate of CCSN-rate and star formation rate 
density in ellipticals as a function of mass. 
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Figure 4. 36 z < 0.2 elliptical Galaxies found to host rpeax < 21.8 CCSNe from the SDSS-II Supernovae Survey, as shown in 30" x 30" coadded images from 
the IAC Stripe 82 Legacy Survey. Crosses (question marks) denote the locations of CCSNe with Pj, < 0.05 (0.05 < Pr, < 0.5). Black (red) symbols denote 
those with (without) spectra. Numbers indicate SN catalogue IDs from the SDSS-II Supernova Survey. 


The black arrow shows the approximate upper limit to counts 
for 10?? — 1095 M,. This is the appropriate one-sided lc error on 
bins with zero counts ( ), normalised for the SN survey 
coverage, TX V, as well as the mean detection efficiency and redshift 
for the aforementioned mass range. 

Additionally there are 2 main sources of systematic uncer- 
tainty on this result. The first relates to uncertainties in SN classi- 
fications within the SN Sample: The CCSN-rate density described 
above is formed from the SN sample excluding those with Pj, > 
0.05. Our positive systematic uncertainty equates to the increase to 
densities found when including 0.05 < Pj, < 0.5 CCSNe. 


A second systematic relates to the dust extinction assumed 
within the host galaxies, or more specifically, the dust screen in the 
region of the SNe. Equations 5 and 6 rely on an estimate of this 
extinction level. From the full CCSN sample of 
( ), the mean r-band extinction in the SN regions was esti- 
mated at 0.5 mag. We therefore also assume a mean extinction of 
0.5 mag for elliptical hosts in the present work. However, as the 
sample in ( ) consisted mostly of star-forming 
galaxies, which can be abundant in dusty star-forming regions, the 
level of dust attenuation in elliptical hosts may be lower. Indeed. the 
median Galactic-extinction corrected absolute r-band magnitude of 
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Figure 5. The core-collapse supernova rate density (left-hand vertical axis) and the implied star formation rate density (right-hand vertical axis), as a function 
of galaxy stellar mass. The red solid line depicts the contribution to these densities from elliptical galaxies. The solid filled red region depicts the associated 
statistical uncertainty, which is the quadrature sum of MC + Poisson + Cosmic Variance errors. Hatched red regions depict systematic errors, from the 
treatment of 0.05 < Pr, < 0.5 objects (positive error) and from the dust extinction assumed for the elliptical hosts (negative error). The blue (black) line and 
filled region shows the CCSN-rate and SF-rate as a function of mass for blue star-forming disk galaxies (all galaxies). The black arrow shows the upper-limit 


at 10°? — 10?5M, at which number counts are consistent with zero (see text). 


Py, < 0.05 CCSNe at peak light is —17.70 mag in ellipticals, com- 
pared with —17.39 mag in blue disks. This implies a median r-band 
extinction of ~ 0.2 mag in ellipticals. For comparison, we find a dif- 
ference of 0.16 mag in median A, estimates when moving from blue 
disk Type Ia SN hosts (median of A, = 0.38 mag) to elliptical Type 
Ia SN hosts (median of A, = 0.22 mag), where V-band extinction 
values were estimated by Sako et al. (2018) from a MLCS2k2 SN 
light curve fitting technique (Jha et al. 2007). Despite these num- 
bers we use a more conservative negative systematic uncertainty on 
the elliptical SFRD, which represents the decrease to densities as- 
suming zero extinction in the SN regions. Systematic uncertainties, 
(labelled ‘sys’) are shown as the hatched regions in Figure 5. 

The SFRD in ellipticals is constrained above zero for a mass 
range 10? - 10!!? Mo. The peak contribution most likely comes 
between 10/96 - 10/93 Mo, where ellipticals contribute 1.6 + 0.8 
(stat) *03 (sys) x10-> CCSNe yr“! Mpe™ dex"! h3,, corresponding 
to an SFRD of 1.3 + 0.6 (stat) *02 (sys) x10? Mo yr! Mpc? dex"! 
hi, at these masses. All the results of the present work are derived 
using a 737 cosmology (h = 0.7, Qm = 0.3, Qa = 0.7). 

There is a narrow range of masses which contribute non- 
negligibly to the SFRD in ellipticals, with ~50% of the contribution 
coming for 10/94 - 10!°8 Mo, and ~90% of the contribution found 
between 10/9? - 10!!? Mo. The integrated density above 10^? Mg, 
corresponding to what is effectively the total elliptical population, 
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is found to be 1.3 + 0.3 (stat) 1: (sys) x10? CCSNe yr! Mpc? 
h3,, equating to 1.0 + 0.3 (stat) *07 (stat) x10? Mo yr! Mpe™? hå. 

For a comparison with the elliptical SFRD we show the result 
from the sample of CCSNe found in star-forming ‘blue disk’ galax- 
ies (shown in blue in Figure 5 and in successive figures), which we 
define as those z < 0.2 galaxies which have (g — r)o < 0.65 mag, 
fav < 0.5, and were found to be non-ellipticals using our classifi- 
cation procedure. These galaxies are expected to be the dominant 
constituents of the star-forming main sequence. 

The SFRD in blue disks as integrated above 10!°° Mọ is a 
factor 2.1 + 0.7 (stat) *13 (sys) times that in ellipticals. The SFRD 
of ellipticals most likely surpasses that in blue star-forming disks 
above 10!°° Mo, due mainly to the domination of the total galaxy 
number density from ellipticals at the highest galaxy masses. 

We also show, in black, the total SFRD as determined from 
the SDSS-II SN Sample. Note that the red and blue series do not 
sum to give the total SFRD: There will be additional contributions 
from lenticular galaxies, irregular galaxies and dusty star-forming 
galaxies, which fall in to neither the elliptical or blue disk groups. 
Our results imply that elliptical galaxies contribute 11.2 + 3.1 (stat) 
I (sys) 96 of the total star-formation budget at present epochs. 
This rises to 20.2+6.0 (stat) *37 (sys) % for masses above 10'°° Mo. 
Both values are consistent with that of Kaviraj (20142), who find 
that 1446 pertains to early-type galaxies. We conclude that ellipti- 
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Figure 6. The galaxy stellar mass function. Solid lines and filled regions 
show z < 0.2 galaxy number densities, and their 1 c MC+Poisson+Cosmic 
Variance errors, respectively, derived from a 1/Vmax method using r < 20 
mag galaxies. In red: elliptical galaxies. In blue: ‘blue disk’ galaxies (see 
text). In grey: galaxies fitting neither classification. In black: all galaxies. 


cals contribute a non-negligible fraction of the total star-formation 
budget to the present-day cosmic volume. 


5 THE SPECIFIC STAR FORMATION RATES OF 
ELLIPTICALS 


Equation 2 demonstrates that the mean SSFR of ellipticals can be 
derived using our result for the SFRD if we also have a calculation 
of the GSMF of ellipticals. To calculate the required GSMF, we use 
a 1/Vmax method applied to the sample of 27940 Bayesian classified 
ellipticals with z « 0.2 and r « 20 mag. The maximum luminosity 
distance, Dt max, within which each galaxy would remain brighter 
than our flux limit of r = 20 mag is given by 


5 log(Dz max) =5 log(Dz obs) +20- lauro + dK , (8) 


where dK is the k-correction for our galaxy at the observed red- 
shift, minus that at the maximum redshift, and where the required 
maximum redshift is iteratively inferred starting from the value of 
Dimax With dK = 0. Finally Dr may values are clipped to lie at or 
below the inferred luminosity distance at z = 0.2. D, max leads to 
1/Vmax, the values for which act as weights on galaxy number den- 
sities in log(M/M.o) bins of width 0.2 dex, yielding volume-limited 
results. We calculate the GSMF of blue disks within our redshift 
and magnitude limits (a total of 66288 blue disk galaxies), and that 
of all z < 0.2, r < 20 mag galaxies (a total of 113239 galaxies), 
using the same approach. 

Corresponding results are shown in Figure 6. Similar to the 
masses at which the peak of CCSN production is observed, the 
number density of z < 0.2 ellipticals peaks between galaxy stellar 
masses of 10!°4 and 10/93 Mo, at 2.11 + 0.03 x 10? dex"! Mpc? 


hy Ellipticals dominate the total z < 0.2 galaxy number density 
for masses > 10/93 Mo. 

Figure 7 shows (with the red solid line and filled regions) the 
mean SSFR of ellipticals as a function of mass, derived from the 
results for the volumetric SFRD and the GSMF. Using an inverse- 
square-error weighting on each mass bin, we find that the mean 
SSFR of ellipticals as a function of mass is well represented by the 
following regression line for masses > 101° Mo: 


log(S(M)/yr) = — (0.80 + 0.59) log(M/10!°°M,) —10.83 + 0.18 . 
(9) 


The uncertainty on the gradient for this regression line im- 
plies there is an 1846 chance that S is constant with mass for ellip- 
ticals. In relation to this uncertainty, we note that although the esti- 
mate of mean SSFR in ellipticals approaches that of main sequence 
galaxies at log(M/Mo) ~ 10.0, the uncertainty on the result for 
individual bins is such that S at log(M/Mg) = 10.0 could be as 
low as 107!!? yr-! within 1c errors, in line with the mean for the 
elliptical population. 

Given uncertainties on the SSFR for ‘blue disks’, which 
effectively define the star-forming main sequence, it is unclear 
whether the result shows a step function at M ~ 10°°Mo and is 
then flat for higher masses, or whether the SSFR instead exhibits 
a more gradual negative slope with mass. A slope of some form is 
likely, given that more massive star-forming galaxies formed earlier 
and a longer time has passed since the peak of their star formation 
(see, e.g. Gallazzi et al. 2005). 

There may however be reason to expect a step function in 
SSFR. For instance, McGaugh et al. (2017) find a distinct star 
forming main sequence for 8.0 < log(M/Mg) < 10.0 compared 
to higher masses, with results consistent with constant S in this 
mass regime. Type II SN-rate models of Graur et al. (2015), which 
build on the work of Li et al. (2011), imply a step function in SS vs 
mass centred on log( M/M;) ~ 9.5. Furthermore, it was found by 
Sedgwick et al. (2019b) that assuming a constant SSFR for masses 
< 10??M, yields the best consistency between the GSMF from 
CCSNe and that from a 1/Vmax method. Future larger samples of 
SNe utilised with these SN-based methods will allow us to reach a 
conclusion for the presence of a step function in S versus mass for 
main-sequence galaxies (e.g. Ivezić et al. 2019). 

The value of S(M) for the total galaxy population could 
be argued to follow a Schechter parameterisation (Schechter 
1976). Fitting for M > 10!°°Mo, we find best-fit parameters of 
[log(M*/Mo), S*/10! yr, a] = 
[10.90 + 0.08, 7.94 + 2.54, —1.41 + 0.03]. 

Similar to Figure 5, the black arrow shows the 1c upper-limit 
on S for 109? — 10?5M,, ellipticals. This value is derived from the 
one-sided 1c error on bins shown in Figure 5, and adopting the av- 
erage galaxy number density and mass in the aforementioned mass 
range for the conversion in Equation 2. This upper limit shows that, 
because of the form of Equation 2, a fractionally large Poisson er- 
ror on host galaxy counts leads to a particularly large error on the 
SSFR at low elliptical masses, spanning several dex. To assess S 
in ellipticals using CCSNe for M < 10?^M, would require future 
data from SN surveys of longer time-span and/or larger sky cover- 
age (e.g. Ivezić et al. 2019) than those currently available. 

Points in Figure 7 show independent SSFRs for individ- 
ual galaxies derived from UV/Optical SED fitting applied to the 
GSWLC-2 sample. 7714 of 27940 (8502 of 49641) ellipticals (blue 
disks) are matched with GSWLC. 12 of 27 (83 of 174) elliptical 
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Figure 7. Left: The specific star formation rate as a function of galaxy stellar mass (via multiple means). The red solid line and filled region denote the mean 
and 1 c uncertainty (MC + Poisson + Cosmic Variance) on the specific star formation rate in elliptical galaxies at z < 0.2. Hatched-regions depict additional 


systematic uncertainties described in the caption of Figure 5. Small-red points depict the SED-derived specific star formation rates of those individual elliptical 
galaxies cross-matched with GSWLC-X2 data. Large-red points denote cross-matched elliptical CCSN hosts, with error bars showing 1 c errors from SED- 
fitting. The blue series depict the same but for blue star forming galaxies. The black series shows the mean specific star formation rate as a function of mass 
for all galaxies. The black arrow shows the upper-limit at 10°? — 10°°Mo at which SN counts are consistent with zero (see text). Right: The same results as 
the left-hand panel but derived using the SFRD and GSMF both integrated over galaxy masses > 1019? Mo. 


CCSN hosts (blue disk CCSN hosts) are matched with GSWLC-2, 
where we here define CCSN hosts as those with Pra < 0.05. Hosts 
are shown as larger points, with their 1o errors on Ssgp. 

We find that there are no significant differences in the opti- 
cal colour distributions of the GSWLC-matched samples compared 
with the full samples of ellipticals and blue disks. We interpret this 
as a lack of evidence for a bias towards the selection of higher Ssgp 
galaxies in the GSWLC-matched sample. As such, a comparison of 
the CCSN-derived and SED-derived results is constructive. 

Ssgp values are qualitatively consistent with S(M) derived 
from CCSNe, with 10 of the 12 GSWLC-matched elliptical CCSN 
hosts having a measurement of Ssgp within the 1o uncertain- 
ties on the SN-derived S(M), and with the running mean of all 
GSWLC-matched ellipticals lying within a lo separation from the 
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SN-derived mean as a function of mass, as shown in Figure 8. 
Interestingly, the running mean value of Ssgp for the main se- 
quence shows no evidence for the aforementioned step function 
about M = 10°°Mo. Furthermore, the gradient of Ssgp with mass 
is comparable for both ellipticals and blue disks, which may be in 
contrast to our findings using CCSNe. 

Integrating the volumetric SFRD and the GSMF for M > 
10!°°Mo (separately) leads to the mean SSFR over this galaxy mass 
range. Mean values are indicated in the right-hand panel of Fig- 
ure 7. From a CCSN technique, we obtain a value of S = 9.2 + 2.4 
(stat) 1d (sys) x107? yr“! for M > 10!°°M, elliptical galaxies. 

Also shown is the counterpart value of Sskp. taken as the 
average SED estimate (in linear space) for all GSWLC-matched el- 
lipticals with M > 10!9?M,. We find Ssgp = 5.7 x 1072 yr! for 
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Figure 8. The specific star formation rate as a function of galaxy stellar 
mass (via multiple means). The red and blue solid lines and surrounding 
filled regions denote the CCSN-derived results and their uncertainties for 
elliptical galaxies and blue disk galaxies, respectively, as described in Fig- 
ure 7. The red and blue dotted lines show the running mean of specific star 
formation rates from SED fitting (derived by Salim et al. 2018), for ellipti- 
cal and blue disks, respectively, averaged over 0.1 dex in mass. Surrounding 
filled regions show the corresponding standard error of these means. 


all GSWLC-matched ellipticals, and Ss =4.1x10-? yr! for the 
sub-sample of GSWLC-matched elliptical CCSN hosts. Assuming 
an average dust extinction in hosts of A, = 0.5 mag, these results 
are in 1.40 tension, and 2.1c tension, respectively, with the value 
of S measured using a CCSNe. Assuming instead zero dust extinc- 
tion in elliptical hosts, we find respective discrepancies of 0.50 and 
l.lo. 

The mean SSFR result for M > 10!°°Mz blue disks, derived 
using a CCSN-based method comes to S = 7.8 + 1.3 x10-!! yr-!. 
This is in < lo tension with SED-based measurements, found to 
be Ssep -9.1x 107! yr! for all GSWLC-matched blue disks, and 
Ssgp = 7.9x 107! yr! for all matched blue disk CCSN hosts. This 
consistency reassures our faith in the elliptical result, given that the 
mean SSFR of the main sequence is well-defined from SED-based 
measurements (and other methods) within the literature (see, e.g. 
Noeske et al. 2007; Speagle et al. 2014). The CCSN-derived mea- 
surement of S for M > 10!°°M. main sequence galaxies implies 
that ellipticals have a mean SSFR which is 11.8+3.7 (stat) +33 (sys) 
% of that on the star-forming main sequence. 


6 THE MEDIAN SPECTRUM OF ELLIPTICAL CCSN 
HOSTS 


Note that the SED result averaged over elliptical CCSN hosts is 
comparable to that averaged over all ellipticals. This implies that 
whilst CCSNe hosts are indeed probabilistically determined by 
their star formation rates, the selection function is not so extreme 
that we are simply tracing the very highest star formation rate el- 
lipticals. This increases confidence that the measurement of S via 
CCSNe is well-representative of the mean star formation level over 
the total elliptical population. 


Of course, SED-based measurements of star formation rates 
for ellipticals are subject to uncertainty; which is indeed a major 
motivation for the present work. We therefore turn to an analysis 
of the 25 elliptical CCSN host galaxy's coadded spectra, to test 
whether these hosts exhibit spectral properties typical of the aver- 
age elliptical. 

Although SED measurements are derived in part from this 
spectral information, we can observe a more detailed picture of the 
star formation properties from the direct spectra, and can look at 
information free of fitting dependencies. 

We first compute the median spectrum of the Pra < 0.05 ellip- 
tical CCSN host galaxies for which spectra are available, i.e. those 
25 galaxies labelled with black crosses in Figure 4. To do this, each 
galaxy's rest-frame spectrum is normalised to equal unity when in- 
tegrated over wavelengths of 4000À and 8000À, before the sam- 
ple's median normalised flux is taken as a function of wavelength 
axis at intervals of 0.1À. The 1 c error on the median flux is calcu- 
lated from the 16" and 84" percentiles as a function of wavelength. 
This uncertainty is shown in each panel except the top main panel, 
where it is omitted for clarity. 

We similarly compute the median spectrum of a control sam- 
ple of non-hosting ellipticals matched in galaxy stellar mass and 
redshift to our CCSN hosts. For each host, we find the separation 
with non-hosts in a Cartesian mass-redshift space, where separa- 
tions along each axis are normalised by the median separation of 
random pairs, such that mass and redshift carry comparable weight 
in the matching procedure. We find the nearest 5 matches for each 
host elliptical, yielding a control sample of 125 galaxies. 

As the Balmer lines are notable signatures star-formation 
properties, we give particular focus to the strengths of Ha and Hf. 
We also investigate the ratios of these lines with the respectively 
nearby [NII]46583 and [OIII]45007 lines, in order to disentangle 
flux contributions to the Balmer lines from star formation with 
those from nuclear activity (Kauffmann et al. 2003). 

As shown in the bottom panel of Figure 9, we find no sta- 
tistically significant difference in Ha or Hf line strength between 
the SN host sample and control sample. We also find no statis- 
tical difference between the median strengths of [OIII]45007 or 
[NII]46583. In fact, not one optical emission line is found to differ 
in its median strength over the 2 samples, within the errors shown. 

We also see no difference in the median line ratios 
[OIII]45007/Hg or [NIT]46583/Ho, implying a similar level of both 
star formation and nuclear activity for the host sample and control 
sample. The equivalent width of these 4 emission lines is also con- 
sistent across the 2 samples. 

We finally use these line ratios to consider the positions of 
the 25 elliptical CCSN hosts on the BPT diagram (Baldwin et al. 
1981), and find all 25 of these galaxies lie comfortably on the 
AGN/Passive side of the Kauffmann et al. (2003) demarcation. This 
implies the line ratios are incompatible with high-levels of star- 
formation and instead the majority contribution to the aforemen- 
tioned lines is coming from old stellar populations and/or AGN 
(which would be Type-II AGN in the majority of cases, given the 
lack of broad lines seen in Figure 9). 

Each of these results implies once more that the elliptical 
CCSN hosts have typical levels of star formation when compared 
to the average (or classical interpretation of an) elliptical, and as 
such we have derived an accurate measure of the mean SSFR of all 
elliptical galaxies, directly from CCSNe. The consistency between 
the spectral properties of the CCSN hosts with the remaining ellip- 
tical population may also imply that the star formation in ellipticals 
is not particularly bursty in nature. 
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Figure 9. Top panel: The median spectrum of elliptical CCSN hosts (red) vs that of a control sample of (non-hosting) elliptical galaxies (see text). Inset plots 
show zoomed in wavelength ranges about the H6 and Ha lines. Upper and lower bounds on the median normalised flux are signified by filled regions and 
correspond to the 84"" and 16" percentiles, respectively. These are omitted in the top main panel, for clarity. Middle panel: The residual spectrum (host - 
control). Bottom panel: The statistical significance of the difference between the median spectra as a function of wavelength. 


7 SUMMARY 


We have classified all z « 0.2, r « 20 mag Stripe 82 galaxies using 
a Bayesian method, trained on the colour and morphology proper- 
ties of manually-classified samples of ellipticals and non-ellipticals 
from Galaxy Zoo 1 (Figure 1). Such a concept could certainly be 
extended to classify other galaxy types, including spiral, lenticular, 
and irregular galaxies, efficiently over large surveys. In doing so, 
the connections between Hubble Type and galaxy evolution could 
be probed with never-before-seen statistical vigour. 

We have isolated a sample of 421 r « 21.8 mag, z « 0.2 likely 
CCSNe, originating from the SDSS-II Supernova Survey. From 
their light curves, these SNe each had a Bayesian confidence (Pcc) 
> 50% of being a CCSN. 360 (~ 86%) of these had a Bayesian con- 
fidence > 95%, making this a reliable sample of CCSNe. 24-196 of 
these CCSNe are likely to be of Type Ib/c, in line with expectations 
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for the ratio of those of Type Ib/c to Type II. As it is far more com- 
mon to misclassify Type Ia SNe as Type Ic SNe (lack of light curve 
plateau and similar spectral properties) rather than as Type II SNe, 
this reasonable ratio of Type Ib/c SNe to Type II SNe increases our 
faith in the SN classification procedure. In Sedgwick et al. (2019b), 
these SNe were carefully matched to their host galaxies, meaning 
we could here isolate those CCSNe residing in ellipticals. 36 CC- 
SNe (27 with Pcc > 0.95) were found to occur in such galaxies 
(Figure 4). 

The z < 0.2 volumetric star formation rate density (SFRD; 
Figure 5) in elliptical galaxies was calculated from the CCSN 
rate as a function of host galaxy stellar mass, assuming a ratio 
between core-collapse events and star formation of R = 1077? 
M5. The elliptical SFRD result was constrained above zero for 
9.8 < log(M/M,) < 11.2, with the peak star formation contribu- 
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tion from ellipticals coming at 10.4 < log(M/M.,) < 10.8. Ellipti- 
cals were found to dominate the present-day star formation budget 
for log(MU/Ms) > 10.9. 

Ellipticals were found to contribute 11.2 + 3.1 (stat) *29 (sys) 
% of the total star-formation budget. This rises to 20.2 + 6.0 (stat) 
13 (sys) % for log(M/Ms) > 10.0. These results are consistent 
with recent studies of star formation in ellipticals using galaxy 
emission, such as the UV-optical study of (Kaviraj 2014a). As 
CCSN statistics cannot be misconstrued as anything but signatures 
of recent star formation, we have definitively demonstrated that e/- 
lipticals contribute a non-negligible level of star formation to the 
present-day cosmic budget. 

The mean specific star formation rate (SSFR; Figure 7) of 
ellipticals was calculated as a function of mass, via the com- 
bined information of the previous SFRD result and the galaxy 
stellar mass function (Figure 6). The mean elliptical SSFR most 
likely exhibits a slope, with a regression line well-representative for 
galaxy masses > 10!°°Mo, such that log(S(M)/yr) = — (0.80 + 
0.59) log(M/10!°°M,) — 10.83 + 0.18 . However, these errors 
show there is still an 18% chance the SSFR is constant with mass 
in ellipticals. The mean SSFR for all log(M/Məo) > 10.0 ellipticals 
was found to be S = 9.2 + 2.4 (stat) +27 (sys) x107!? yr7!. This is 
11.8+3.7 (stat) ET (sys) % of the average SSFR on the star-forming 
main sequence. 

An independent mean SSFR in ellipticals derived from 
GSWLC-2 SED fitting is found to be Ssep = 57x 10°? yr, 
which is moderately consistent with our SN-derived result (1.40 
tension). Assuming zero dust extinction in elliptical CCSN hosts 
this tension drops to only 0.5. This indicates the SED-derived re- 
sults put forward as upper limits by (Salim et al. 2018) are likely 
close to the true levels of star formation for these objects (See Fig- 
ure 8). 

We finally computed the median optical spectrum of the 27 
Pia < 0.05 CCSN hosts and compared this with that of a control 
sample of non-hosting ellipticals matched in galaxy stellar mass 
and redshift (Figure 9). We find no statistically significant differ- 
ence in the median strengths of emission lines for these 2 sam- 
ples, including various commonly associated with star formation 
(Ha, HB, [OIII]45007 and [NII]A46583). This result implies ellipti- 
cal CCSN hosts have typical levels of star formation compared to 
the average elliptical. This in turn signifies that our CCSN-derived 
results are an accurate representation for the total elliptical popu- 
lation at current epochs. These results are consistent with the hier- 
archical evolution widely accepted under a A-CDM paradigm, in 
which low-level star formation is expected to continue in ellipticals 
out to current epochs, likely due to the influence of galaxy mergers 
and/or a slow-rate of gas cooling. 

Note that out of all red sequence galaxies, CCSNe were only 
likely to occur in massive red ellipticals within the time-frame of 
the SDSS-II Supernova Survey, and so cannot reveal information 
on the star formation properties of lower mass, environmentally 
quenched galaxies. These may indeed have zero levels of star for- 
mation. A relevant investigation may become possible with access 
to future high-cadence surveys with increased time-span, area and 
magnitude depth, such as LSST. These surveys would greatly in- 
crease the size of CCSN-selected galaxy samples, and would also 
allow for more precise measurements of the SFRD and SSFR in 
massive ellipticals. A study such as this with LSST may offer up 
enough CCSN statistics to allow volume-averaged star formation 
properties to be estimated as a function of redshift, which would 
shed light on the average star formation history of ellipticals. 


8 ACKNOWLEDGEMENTS 


TMS acknowledges support from an STFC DTP studentship, 
jointly supported by the Faculty of Engineering and Technology at 
LJMU. SK acknowledges support from the STFC [ST/S00615X/1] 
and a Senior Research Fellowship from Worcester College Oxford. 

This publication uses data generated via the Zooniverse.org 
platform, development of which is funded by generous support, in- 
cluding a Global Impact Award from Google, and by a grant from 
the Alfred P. Sloan Foundation. The authors acknowledge the role 
of citizen science in the production of Galaxy Zoo data. 

The construction of GSWLC was funded through NASA 
award NNX12AE06G. 

Funding for the SDSS-II, SDSS-III and SDSS-IV has been 
provided by the Alfred P. Sloan Foundation, the National Sci- 
ence Foundation, the U.S. Department of Energy Office of Science, 
the National Aeronautics and Space Administration, the Japanese 
Monbukagakusho, the Max Planck Society, the Higher Education 
Funding Council for England, and the Participating Institutions. 
SDSS-IV acknowledges support and resources from the Center for 
High-Performance Computing at the University of Utah. The SDSS 
web site is www.sdss.org. 

SDSS-II, SDSS-III and SDSS-IV are managed by the Astro- 
physical Research Consortium for the Participating Institutions of 
the SDSS Collaboration including the American Museum of Natu- 
ral History, University of Arizona, Astrophysical Institute Potsdam, 
University of Basel, the Brazilian Participation Group, Brookhaven 
National Laboratory, University of Cambridge, the Carnegie Insti- 
tution for Science, Carnegie Mellon University, Case Western Re- 
serve University, University of Chicago, the Chilean Participation 
Group, University of Colorado Boulder, Drexel University, Fermi- 
lab, University of Florida, the French Participation Group, the Ger- 
man Participation Group, Harvard-Smithsonian Center for Astro- 
physics, Harvard University, the Institute for Advanced Study, In- 
stituto de Astrofísica de Canarias, the Japan Participation Group, 
Johns Hopkins University, the Joint Institute for Nuclear Astro- 
physics, the Kavli Institute for Particle Astrophysics and Cosmol- 
ogy, Kavli Institute for the Physics and Mathematics of the Uni- 
verse (IPMU) / University of Tokyo, the Korean Scientist Group, 
the Chinese Academy of Sciences (LAMOST), Lawrence Berke- 
ley National Laboratory, Leibniz Institut für Astrophysik Potsdam 
(AIP), Los Alamos National Laboratory, Max-Planck-Institut für 
Astronomie (MPIA Heidelberg), Max-Planck-Institut für Astro- 
physik (MPA Garching), Max-Planck-Institut für Extraterrestrische 
Physik (MPE), the Michigan State/Notre Dame/JINA Participa- 
tion Group, Universidad Nacional Autónoma de México, National 
Astronomical Observatories of China, New Mexico State Univer- 
sity, New York University, University of Notre Dame, Observatário 
Nacional / MCTI, Ohio State University, University of Oxford, 
Pennsylvania State University, University of Pittsburgh, University 
of Portsmouth, Princeton University, Shanghai Astronomical Ob- 
servatory, the Spanish Participation Group, University of Tokyo, 
the United States Naval Observatory, University of Utah, Vander- 
bilt University, United Kingdom Participation Group, University of 
Virginia, the University of Washington, University of Wisconsin, 
Vanderbilt University, and Yale University. 


9 DATA AVAILABILITY 


The data underlying this article were accessed from: the Sloan 
Digital Sky Survey at skyserver.sdss.org (dr14.PhotoPrimary, 
dr14.SpecObj, dr14.Photoz, dr7.PhotoObjAll) and at data.sdss. 


MNRAS 000, 1—15 (2021) 
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